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Abstract 



The generalized quark-antiquark potential of A^ = 4 supersymmetric Yang-Mills 
(D ' theory on S'^ x M calculates the potential between a pair of heavy charged particles 

separated by an arbitrary angle on §^ and also an angle in flavor space. It can be 
calculated by a Wilson loop following a prescribed path and couplings, or after 
>■ \ a conformal transformation, by a cusped Wilson loop in fiat space, hence also 

generalizing the usual concept of the cusp anomalous dimension. In ArfS's x S^ 

\Q \ this is calculated by an infinite open string. I present here an open spin-chain 

model which calculates the spectrum of excitations of such open strings. In the 

w ■ dual gauge theory these are cusped Wilson loops with extra operator insertions 

^ ' at the cusp. The boundaries of the spin-chain introduce a non-trivial reflection 

phase and break the bulk symmetry down to a single copy of psu(2|2). The 
dependence on the two angles is captured by the two embeddings of this algebra 

^ ! into psu(2|2)^, i.e., by a global rotation. The exact answer to this problem is 

j^ \ conjectured to be given by solutions to a set of twisted boundary thermodynamic 

Bethe ansatz integral equations. In particular the generalized quark-antiquark 
potential or cusp anomalous dimension is recovered by calculating the ground 
state energy of the minimal length spin-chain, with no sites. It gets contributions 
only from virtual particles reflecting off the boundaries. I reproduce from this 
calculation some known weak coupling perturtbative results. 




Figure 1: A cusped Wilson loop (a) in A/" = 4 SYM on M^ (and by analytic continu- 
ation also on M.^'^) is related to a pair of antiparallel lines (b) on §^ x R. In the dual 
string theory this is calculated by a string world-sheet in AdS^ x S^ ending along 
the two lines on the boundary (c). 



1 Introduction 



Wilson loops are some of the most important observables in nonabelian gauge theories. 
Among their many features, they capture the potential between heavy probe particles, hence 
can serve as an order parameter for confinement [1] . They calculate a big part of the effect of 
high-energy scattering of charged particles [21 13], and in the case of A/" = 4 supersymmetric 
Yang-Mills theory (SYM) are conjectured to calculate scattering amplitudes exactly |11|5]. 

In [6] a family of Wilson loop operators in A/" = 4 SYM were presented and studied. Two 
useful points of view on these operators are (see Figure [1]) 

• A cusp of angle vr — in W^, i.e., two rays meeting at a point. 

• Two lines along the time direction on §^ x M, separated by an angle vr — (/) on the sphere. 

These two configurations are related to each other by a conformal transformation, so are 
essentially equivalent. The Maldacena- Wilson loops of A/" = 4 SYM also include a coupling 
to a real scalar field and it is natural to allow the two rays/lines to couple to two different 
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scalars, say $^ and $^ cos 6* + $^ sin 9. Thus giving a two parameter family of observables, 
with the opening angle vr — </> (such that </> = is the straight line) and 9. 

These geometries calculate the physical quantities mentioned above. The pair of antipar- 
allel lines gives the potential V{\, 0, 9) between two heavy charged particles propagating 
over a large time T on S^ 

(Wynnes) = e-^^(^'^'^) . (1.1) 

It is possible to recover the potential in fiat space by taking residue of ^ as — )■ vr (see [6], [7] 
for a detailed discussion). Cusped Wilson loops have the structure 

(W^cusp) = e-^--(^'*''')i°^T, (1.2) 

with the same function rcusp(A, 0, 6') = y(A,0, 6') and A and e are IR and UV cutoffs, 
respectively. Scattering amplitudes are related to cusped Wilson loops in Minkowski space, 
so one should take (p = iip to be imaginary. In particular in the limit of (/? — ?■ oo, leads to 
7cusp — the universal cusp anomalous dimension 

r_p(A,z(^,^)^^^;^. (1.3) 

As expalined in detail in [8], in the limit of small and vanishing 9, rcusp(A,0, 0) ~ 
—B{X)(f)^ calculates also the radiation of a particle moving along an arbitrary smooth path. 
An exact expression is given there for B including 1/A^ corrections to all orders in the gauge 
coupling A. A modification of this expression applies also in expanding Fcusp around the BPS 
configurations = ±6* [9j. In this manuscript an integrable spin-chain model is presented 
which is conjectured to calculate the full function Fcusp in the planar approximation. 

There is another application of this quantity, where it was originally defined [TOl [TT] . 
Wilson loops satisfy a set of equations known as the loop equations [12]. These equations 
are nontrivial only for intersecting loops, where cusps appear. The loop equations have been 
solved in zero dimensional matrix models and in two dimensional Yang-Mills. To define 
them properly in four dimensional theories requires understanding cusped Wilson loops and 
how to renormalize their divergences. Logarithmic divergences of Wilson loops (and in the 
case of the Maldacena- Wilson loops in A^ = 4 SYM, all the divergences) come from cusps, 
or from insertions of local operators. As explained in the following, integrability supplies the 
answer to this question in A/^ = 4 SYM, where all the logarithmic divergences arising from 
either insertions or cusps are governed by the same spin-chain model (to be precise, a set of 
twisted boundary thermodynamic Bethe ansatz equations (TBTBA, BTBA, or TBA in the 
follwoing) ) . 

Integrability in this calculation is most easily seen in the string theory dual, where at 
leading order the Wilson loop is calculated by a classical string solution [131 [IS [ISl CSl E] 
and of course the string sigma-model in AdS^ x §^ is integrable [17]. At the one-loop 
level all the fluctuation operators about these particular classical solutions turn out to be 
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Figure 2: A generalization of Figure [T] allows an arbitrary adjoint valued local oper- 
ator to be inserted at the apex of the cusp (a) in M^. After the conformal transfor- 
mation this is a pair of antiparallel lines (6) on S^ x M, seemingly like in Figure [1^, 
but in fact the details of the operator O are represented by nontrivial boundary 
condition \i\}) at past and future infinity. The dual string solution in AdS^^ x S'^ 
still ends along the two lines, but is in an excited state (in the figure the sphere is 
suppressed) . 



versions of a known integrable system, the Lame differential operator [181 [IHl 1201 [6]. The 
calculations of the same quantity on the gauge theory side were based on brute force Feynman 
diagrammatics. 

In order to see the integrability of these operators also in the gauge theory, it proves 
useful to consider a further generalization of the set of observables. Along a Wilson loop in 
the fundamental representation one can insert local operators in the adjoint representation 
of the gauge group while retaining gauge invariance. We shall insert such an operator at 
the apex of the cusp (Figure [2]). In the S"^ x M picture this operator is at past infinity, and 
corresponds to an excitation of the Wilson loopllj 



^The infinite lines are not well defined observables, as they are not invariant under gauge transformations 
at infinity. One should supply extra boundary conditions to construct them, and usually one assumes that 



As mentioned, the expectation value of cusped Wilson loops suffer from logarithmic 
divergences, which can be considered the conformal dimension of the operator, meaning that 
these observable are eigenstates of the dilation operator. The generalized quark-antiquark 
potential V in (11 .ip (and Fcusp in (ll.2p ) is exactly that conformal dimension. In order to 
understand this dimension as the solution of a spin-chain problem it is useful to consider 
the more general situation. 

Note that in order to make this argument it is crucial that geometrically the cusp is 
invariant under dilations and the antiparallel lines under translations (also in addition a 
f/(l) of rotations in the transverse space and an 50(4) flavor symmetry). Therefore one can 
define the problem of finding operators in the quantum theory which are eigenstates under 
this transformation. A similar argument cannot be implemented in a trivial way for Wilson 
loops of other shapes. 

Another way to visualize the situation in the dual string theory, where the Wilson loops 
are described by a string stretching between the two lines at the boundary of AdS^ x S^. 
Thus, along the boundary of the open string we have two different boundary conditions, 
corresponding to the two lines. In world-sheet language we should consider the string with 
the insertion of two operators, these are not regular open string states, rather they are 
boundary changing operators. To be precise, for any pair of boundary conditions there is 
a family of such boundary changing operators exactly corresponding to the cusped Wilson 
loops with extra insertions. The original cusp without an extra operator insertion maps 
to the ground state of the open string, i.e., the boundary changing operator of the lowest 
dimension. 

The case of the insertion of adjoint operators into the straight Wilson line with 6 = (p = 
was studied in [2T]. In the string theory dual these operators are mapped to regular open- 
string vertex operators as the boundaries of the string do not get modified. On the gauge 
theory side the dimension of these insertions can be calculated by solving an open spin-chain 
problem, as is reviewed in the next section and extend to the more general situation. 

If the operators carry macroscopic amount of angular momentum in AdS^ or S^, then 
there is a semiclassical description for them in string theory. For the insertions into the 
straight line this was done in [21], and it can be generalized for an arbitrary cusp using 
the techniques in [22]. These solutions are such that the string approaches the boundary 
along the same two lines as the solutions without the insertions (see Appendix B of [B]), but 
performing extra rotations around the §^ or in AdS^ in the bulk of the world-sheet. 

To be completely explicit, I write down the form of the Wilson loop operator with inser- 



thc lines close onto each other in a smooth way. The excitations considered here have alternative boundary 
conditions. 




Figure 3: Yet another picture of a Wilson loop with two cusps (with possible local 
insertions) connected by arcs. It is related to that in Figure |2]2 by a conformal trans- 
formation mapping the point at infinity to finite distance. If the distance between 
the cusps is d, the expectation value of this Wilson loop is (W) oc 1/(P^^^''^'^\ where 
the logarithmic divergences in f ll.2p are interpreted, as usual as renormalizing the 
classical dimension. 



tion of the simplest adjoint valued local operator O = Z'^ with Z = $^ + i$^ 

_ pg/!J^Mi(x(s))+$i(x(s))ds ^j'gNjg/g°°(icos0yli(x(s))+isin0yl2(x(s))+cose$l(x(s))+sin6»'I«2(x(s))]ds 

x{s) = < (1.4) 

[(scos</), ssin</), 0, 0) s>0, 

By a conformal transformation the point at infinity can be mapped to finite distance (see 
Figure [3]) and then the Wilson loop is a completely kosher gauge invariant operator with 
two cusps with local operators Z"^ inserted at one and Z'^ at the other. 

V{\, 0, 9) was calculated in [6] to second order in perturbation theoryo The result is 
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(lis (e^^^) - C(3) - ^<f> (hi^ (e^^*) + ^) + ^0^) 



In the next section the spin-chain model is presented. Section [3] calculates the leading 
weak coupling contribution to the the cusped Wilson loop from the exchange of a single 
virtual magnon (Liischer-Bajnok-Janik correction). Section H] presents the twisted boundary 
TEA which calculates this quantity at all values of the coupling. Finally Section [5] discusses 
some of the results and possible generalizations. Some technical details are relegated to 
appendices. 

At an advanced state of this project I learnt of [23], which has a great deal of overlap 
with the results reported in this manuscript. 



^Thc three loop answer was derived very recently in [?]■ 



2 Wilson loops and open spin— chains 

In this section the relation between the insertion of adjoint valued operators into a Wilson 
loop and open spin-chains is developed. First the main principles are explained in general 
terms and only later are the precise formulas for the open spin-chain presented. 

In the case of gauge invariant local operators the spin-chain arises [23] by choosing a 
reference ground-state Ti Z^ with Z = <I>^ + z$^ (also known as the BMN vacuum [25]) and 
considering the replacement of some of the Z fields by other fields of the theory: The other 
scalars, the fermions or field-strenghts. In addition one can act with derivative operators. 
The possible insertions are labeled by representations of PS'f/(2|2)^ C P5'f/(2,2|4) which is 
the residual symmetry preserved by the vacuum. These states are viewed as excitation of a 
spin-chain and it is conjectured that this spin-chain is integrable and satisfies a particular 
dispersion relation leading to a solution in principle of the spectral problem of the theory, 
which has passed many stringent tests. The symmetry is reviewed in Appendix |X] and 
the scattering matrix of the fundamental representation of this spin-chain is presented in 
Appendix [Bl 

The same thing can be done with insertions of adjoint valued local operators into Wilson 
loops, as explained some time ago in [21]. That paper studied only the insertions of operators 
made of two complex scalar fields {Z and y = $^ + i$^ forming the SU{2) sector) into a 
straight (or circular) Wilson loop. Here the construction is generalized to the full set of 
allowed states. Also, the crucial modification of the spin-chain model when replacing a 
straight line with an arbitrary cusp is explained. 

The basic principle of [21] , illustrated in Figure HJ is that the calculation of the Wilson 
loop with local insertions is very similar to that of the dimension of local operators. At one 
loop order in the planar approximation every field in the local insertion interacts directly 
only with its two nearest neighbors. Only the first and last fields interact with the Wilson 
loop itself. This leads naturally to an open spin-chain model, where the bulk hamiltonian 
is the same as the usual one as in the closed spin-chain, but with two boundaries where one 
must specify boundary conditions, i.e., reflection matrices. 

It is important to emphasize that this is based mainly on abstract arguments and on 
very few explicit calculations, which is also the case (but to a lesser extent) with the usual 
spectral problem. It is clear from diagramatics that indeed the bulk spin-chain hamiltonian 
is identical to the closed spin-chain one. The hamiltonian in any case is known explicitly 
only at low orders in perturbation theory, but the scattering matrix is known and therefore 
expected to be the same. The boundary interactions were calculated in [21] only in the SU{2) 
sector and only at one loop, where they are essentially trivial. Below, an exact expression 
for the reflection matrix is proposed, which should capture the boundary interactions at all 
loop order. It is clearly conjectural, based on the assumption of integrability, representation 
theory, boundary crossing symmetry and minimality of the dressing factor (see below). The 
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Figure 4: Sample planar Feynman graph calculations for a Wilson loop with two 
insertion of operators made of five scalar fields at the cusps At leading order there 
are only five free propagators (dashed lines) so the classical conformal dimension 
is five (a). At one loop there are nearest neighbor interactions among the scalar 
fields (like b), which are identical to those of single trace local operators, but also 
interaction between the last scalar and the Wilson loop (c). At this order planar 
graphs connecting the loop to itself are restricted to each of the arcs, they are 
finite and do not modify the conformal dimension. Finally The two arcs interact 
by "wrapping effects", which arise in this example only at six loop order (d). Up 
to this order the dimension of an operator made of the scalars Z and Y does not 
depend on the cusp angles and 9. 



formalism presented does reproduce the results of the perturbative calculation, which lends 
credence to this conjecture. 

There are other manifestations of open spin-chains in A/" = 4 SYM, which apart for the 
Wilson loops is whenever there are D-branes in the dual AdS^ x S^ space, on which the open 
strings can endo 

The most studied example is that of the D3-brane giant graviton, which is a determinant 
operator in the field theory, see [261 EH [28] . A closely related system is that of AdS^ filling 
D7-branes introducing fundamental matter into the theory which can serve as the ends of 
an adjoint valued word, again leading to open spin-chains. 

Another system, that of D5-branes representing 3d defects in the gauge theory (domain 



■^Thc strings describing the Wilson loop are also open, but extend all the way to the boundary of space, 
rather than ending on a D-brane. 



walls with fundamental matter) [29] has a symmetry much closer to that of the Wilson 
loops and has been erroneously thought for a long time not to lead to integrable boundary 
conditions on the spin-chain |30] . This has been corrected recently in [311 [32] , and is very 
useful for the derivation below. 

Before going into the technical details, it is useful to explore some of the features that 
can be extracted from looking at the Feynman diagrams, as in Figure JH 

The modification from the straight line to the cusped loop (or from the circle to the pair 
of cusps as in Figure |3]) is quite easy. The spin-chain has to satisfy two different boundary 
conditions at the two ends. The two rays of the Wilson loop couple to linear combinations of 
the scalars $^ and $^ break the 5*0(6) i?-symmetry to SO (4). The ground state Z"^ breaks 
it further to U{1) rotating $'^ into $^. It also has charge J under the U{1) which rotates <l>^ 
into (j)^. 

The fields Z and Y (also Z and Y) interact with the two segments of the Wilson loop 
in the same way, at the one loop order they have purely reflective boundary conditions, so 
the eigenstates of the one loop hamiltonian of the form '^^^e^^'' Z^~^YZ'^~'' are Neumann 
functions. The fields $^ and $^ (or X and X) interact differently with the two boundaries, 
a linear combination of them (depending on the scalar coupling of the Wilson loop) has to 
vanish at the boundary, so they satisfy partially Dirichlet and partially Neumann conditions 
and the allowed magnon momenta will depend on 6. 

Going to higher order in perturbation theory doesn't change much until order J + 1. 
The range of interaction in the spin-chain is identical to the order in perturbation theory 
and at this order there are graphs which communicate all the way from one side to the 
other (Figure JUi). In analogy with the closed spin-chain these can be called "wrapping 
corrections". For operators made of Z and Y these are the first graphs where the answer 
will depend on the cusp angles and 6. In particular the ground state Z"' has dimension 
A = J up to order A"^ and gets corrections only from wrapping effects (and only for (p ^ 6, 
when the system is not globally BPS [9]). The wrapping corrections are hard to calculate 
directly for large J, but in the spin-chain formalism they are given by Liischer like corrections 
[33| ^^, as discussed in Section |3l Higher order wrapping effects are best captured by the 
twisted boundary TBA equations. 

Note though, that in the case of primary interest, that without the local insertion, J = 0, 
so wrapping corrections contribute at one loop order, double wrapping at two loops and so 
on. In that case the diagramatics are not that hard and indeed these calculations were 
done up to three loop order in [151 ESI [HI [7] (see also [361 EZ])- This system then provides 
an interesting laboratory to study the twisted boundary TBA equations, where the desired 
quantity is the ground state energy, which is the simplest observable in the TBA. This is 
discussed in detail in Section Hj below. 



2.1 Boundary symmetries 

It is now time to start with the detailed analysis of the system. The first question is how the 
symmetry of the theory and of the usual spin-chain model is modified by the presence of the 
boundaries. It is crucial that one can focus on one boundary at a time, since the symmetry 
preserved by each boundary is significantly larger and more restrictive than that preserved 
by both together. 

The p5u(2, 2 1 4) supersymmetry algebra of the theory is written down in Appendix |X] as is 
its breaking to psu(2|2)^ x u(l) by the choice of ground state. Here we describe the breaking 
of the symmetry by the Wilson loop, which furnishes the spin-chain with boundaries. 

A straight Wilson loop in the x^ direction and coupling to the scalar $^ preserves half 
the supercharges, those given by the combinations 

Ql + ^'^'i'p.PabQ''' , S^A - ^o.,!'"' p\bS^, ■ (2.1) 

Here 7^^ are the usual gamma matrices in space-time and pi are those for the S'0(6) flavor 
symmetry. 

These generators close onto an osp(4*|4) algebra. It is possible to include an extra phase 
e*" between the two terms in both sums and the algebra still closes to an isomorphic algebra. 
For a = 7r/2 this is the symmetry of an 't Hooft loop and other values correspond to dyonic 
loop operators. So this is the symmetry preserved by each of the boundaries. 

Consider the Wilson loop with /i = / = 1 and the choice of 5*0(6) gamma matrices 

(2.2) 

Of the two copies of psu(2|2) annihilating the vacuum Z"^ flA.3l) . only one diagonal copy 
survives, once considering the linear combinations (12. ip annihilating the Wilson loop. Those 
are 

Q: + ^{^'r,{^')\Ql , si - i{a'f^{a')\S'^ . (2.3) 

A similar algebra, though with a different real form and embedding into psu(2,2|4), 
namely 05p(4|4, R), is also preserved by defects represented in AdS^ x S^ by D5-branes. As 
is shown in [31j the open spin-chain with such boundary conditions has the same boundary 
symmetry which is the intersection osp(4|4, M) fl psu(2|2)^ = psu(2|2). 

In that case the defect breaks the SO{Q) R symmetry to SO (3) x 5*0(3), while for the 
Wilson loop it is broken to 5*0(5). In that case (as in other realizations of open spin-chains 
in A/" = 4 SYM) there are two choices for the vacuum, depending on which copy of the 
two inequivalent 50(3) is broken by the choice of ground state. In one case the resulting 
boundary is charged under the unbroken 50(3), and therefore carries a representation of 
p5u(2|2). In the other case, it is uncharged and does not carry a boundary degree of freedom. 
In the case of the Wilson loop it is natural to break the 50(5) symmetry to 50(3) and the 
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resulting boundary has no degrees of freedom. This can be seen by exphcit calculation, since 
the straight Wilson loop with an insertion of Z^ is a perfectly good BPS operator and does 
not require other fields "shielding" the local operator from the Wilson loop (as would be 
the case with a vacuum like X'^ with X = $^ + i$^)jj This can also be seen from the dual 
string point of view, where the description of the string ground state [2T] does not require 
breaking the symmetry and the resulting Goldstone bosons of global rotation in AdS^ x S^. 
The spin-chain calculation below can be quite easily generalized to the case of boundary 
degrees of freedom to describe both type of D5-brane open spin-chains. 

2.2 Notations 

Before writing down the spin-chain model which calculates these Wilson loop operators let 
us fix some notations. 

The magnons can be characterized by the spectral parameters x^ . Introducing the pa- 
rameter u, then for a general bound state of Q magnons in the physical domain they satisfy 

x^ g 2g 47r 

For generic values of u G C there are four possible solutions to the above equations with x~^ 
and x~ both outside the unit disk, both inside and with either one outside and the other 
inside. 

The momentum p and energy E of the bound state are 



e'P = ^, E = Q + 2ig(^-—]=.b^+16g^^Jl. (2.5) 



x~ \x+ x~y V 2 

It is useful to introduce the generalized rapidity z defined on a torus with half periods 

wi = 2K , L02 = 2iK' - 2K (2.6) 

where K and K' are complete elliptic integral of the first kind with modulus squared k"^ = 
— 16g^/Q^ and k'"^ = 1 — k^. For real g the first period is real and the second is purely 
imaginary. This torus covers the four copies of the u plane. The spectral parameters, 
momentum and energy are expressed in terms of Jacobi elliptic functions (with modulus k) 
of the rapidity z 



I , , Q fcnz 

X [z] — ' 



Lj f cn z \ V 

— ( ±i) (1 + dn 2;) , p(z) = 2 am z , sin - = sn 2 , E(z) = Qdnz . 

Ag \snz J 2 

(2.7) 



''The X"^ state breaks 50(5) — > 50(4) and would require boundary excitations. But it is not a good 
vacuum, as can be seen by the fact that it does not share any supersymmetry with the Wilson loop. 
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Figure 5: Using the reflection trick a semi-infinite open spin-chain can be replaced 
with a spin-chain on the entire line. The original magnons (a) carry representations 
of p5u(2|2)ij X psu(2|2)i;, and momentum p, which gets reflected to momentum —p. 
In the doubled picture (6) the psu{2\2)ji label is carried by a magnon of momentum 
—p on the right side, which gets scattered off the magnon of momentum p and a 
P5u{2\2)l representation. As usual, the momentum does not get modified by the 
scattering and the magnon with momentum — p continues on the left side, now 
carrying a p5u(2|2)i representation 



For real z both x^ > 1, the momentum is real and the energy positive. Shifting z ^ z = 
2; + Ci;2 is a crossing transformation which replaces x^ — )■ 1/x^ and reverses the signs of both 
the energy and momentum. 

The shift z ^ z = z + 072/2 gives the mirror theory, where x~ < 1 < x~^ and both the 
energy and momentum are purely imaginary. In the mirror theory their roles are reversed 
and one defines real mirror momentum and mirror energy as 



p = iE = k'Q sc z , E = ip 



2i arcsin ( — cd ( 

Aw 



(2. 



2.3 Reflections and open spin— chains 

The Wilson loop breaks the psu(2,2|4) symmetry to osp(4*|4) and the symmetry preserved 
by the ground state of the spin-chain psu(2|2)^ to the diagonal component psu(2|2) with the 
supercharges (12. 3p . To write down the spin-chain system preserving this symmetry, one can 
use the the method of images. This was done for the case of the open spin-chain associated 
to the D5-brane domain walls in |31] . 

The construction of the usual p5u(2|2) scattering matrix is reviewed in Appendix [B| 
following [38]. 

Consider first the right boundary and a semi-infinite spin-chain to its left. Each magnon 
with momentum p carries a representation of psu(2|2)i x psu(2|2)j:j. In the doubled descrip- 
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tion one splits each of the magnons in two, one on the original left ray with momentum p 
and a label of psu(2|2)£ and the image magnon of momentum —p and a label of psu(2|2)j:j 
on the right. 

Each magnon is charged under the three central charges of the extended p5u(2|2): C, 
P and K, related also to the labels p, / and a (see Appendix [B]). These charges are the 
same for the two copies of psu(2|2), and one should specify how they transform under the 
reflection. The reflection of the Psu(2|2)r magnons to the mirror description acts as 

_ _ 2.9 

Q^iQ, S^-iS. 

The extra phase acquired by / is crucial, as for neighboring magnons the two /s have to 
satisfy /i = 6*^^/2. For the magnons on the original (left) segment / is given by the momenta 
of all the magnons to their right, up to the wall. For consistency, the image magnons on 
the right, should have / which is the exponent of minus the momenta of all the magnons to 
their left up to the wall including themselves. The definition of / above is exactly this, after 
the refiection. 

In formulas, enumerating the magnons on the left of the wall 1, ■ ■ ■ , M, they satisfy 
fk = 6*^*^+1/^+1 with Jm = 1- Then continue this pattern to the full real line, where the 
mirror magnons are labeled M + 1, ■ ■ ■ , 2M. From (12. 9p fM+k+i = e.^^^'~'' fhi-k so that 
JM+k = e'P^'-f'+^fM-k+i = fu-k = e'P'^'+'^+^fM+k+i satisfies the correct relation. 

If one considers the two particle scattering matrix, the order of the two entries is reversed 
as are the signs of the momenta. Instead of /2 it will depend on the transformed /i — )■ e*^^ /i = 
^ipi+iP2 j^_ Indeed it is easy to check that the S-matrix f lB.6p . f lB.7p is invariant under the 
transformation (12. 9 p 

§(e«pi+^P2j2, -P2, -Pi) = S(/2,pi,p2) . (2.10) 

The open psu(2|2)^ spin-chain on the left ray is the same as a psu(2|2) spin-chain on 
the full line. The symmetry generators in (12. 3|) are the diagonal components of the original 
p5u(2|2)/^ and the refiection (12.91) of psu(2|2)r. The scattering of the right most magnon 
on the left part off its image, the left most magnon on the right part is then the reflection 
matrix for the open spin-chain. The matrix structure of this reflection matrix is flxed by 
extended psu(2|2) symmetry to be the same as the bulk reflection matrix 

'Xp) = Rf\p)k:{p.-p). (2.11) 



where 

"1fa-p) = Tr' "i (2-12) 



_ ^fAv. -p) 

^cd.yP'! P) — TTT r 



is the scattering matrix with the scalar factor removed. The explicit components of M^^-* (p) 
are written down in Appendix O The scalar part of the reflection matrix is not flxed by 
symmetry, but it is constrained by crossing symmetry, as discussed shortly. 
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Figure 6: The boundary crossing equation equates the two processes in (a). The 
bottom is a single reflection. The top has a bulk scattering (a pair of psu(2|2) 
scatterings, note that not four, since the two copies of psu(2|2) interact only at the 
boundary), and a crossed reflection. In the doubled picture (b) this single scattering 
(bottom) is related to a pair of regular scatterings and one with both particles 
crossed. 



The boundary Yang-Baxter equation is automatically satisfied, since it can be decom- 
posed into several applications of the bulk Yang-Baxter equation for the psu(2|2) chain. It 
was also checked explicitly in |31j . 



2.4 Boundary scalar factor 

Following [39] it is now commonly accepted that the scattering matrix of the bulk spin- 
chain is invariant under crossing symmetry. The same would then be expected also for the 
boundary reflection. In terms of generalized rapidity variables z, the boundary crossing 
unitarity condition states 



z)=S{z,-z)R'^^^\u2-z) 



(2.13) 



where M'-^^ is the reflection matrix in the crossed channel (see Figure [6]). In our case the 
bulk S matrix is made of two psu(2|2) matrices and including all the indices this is 



bb . 



&^ 



-^)&^, 



-z)C^^C^%, 



!,{R)dd 



dd 



(^2 - Z) 



(2.14) 



where C'^'^ is the charge conjugation matrices. 

Deriving the equation for the boundary scalar phase is actually not that hard, using the 
crossing invariance of the bulk scattering matrix. The crossed reflection matrix is the same 
as the doubly crossed psu(2|2) scattering matrix. After including the usual scalar phase it 



13 



is on its own invariant under both the crossings, giving 



(2.15) 
C%.Sg(^ - ^2, -z) = ,r'Z r... S- (-., z) . 



Ro{u2 - Z) ^g <Sdc/^_,, ^\ _ Ro{^2 - Z 



Expressing also the reflection matrix on the left hand side of f l2.14p in terms of the bulk 
scattering matrix and using unitarity gives 

Ro{z)Ro{z-U2) = So{z,-z)So{z-uj2,uj2- z). (2.16) 

This equation can be solved with Ro{z) = So{z, —z). 

This expression, though needs to be properly defined. The reason is that the solution 
of the crossing equation was formulated in terms of the dressing factor a{zi,Z2), which is 
related to 5*0(2:1,22) by 

So{zuZ2r = ^i~V~]{T'- ^^2 (2-17) 

x-^ — X2 1 — 17x^X2 a{Zi,Z2)^ 

and in particular 



^ , , X~^ X + 1/x 1 

So{z,-z) J— T^—, r. 2.18 

^ ^ x-\j x+ + l/x+ a{z,-z) ^ ' 

So in defining S^i^z^ —z) one has to specify how to deal with the square root. 
It is natural therefore to take an ansatz for the boundary phase as 

Ro{z) = ^^ (2.19) 

a{z,-z} 

such that equation (I2.16P is 

(TB{z)aB{z - ^2) = 5*0(2;, -2)5*0(2 - UJ2,(^2 - z)a{z, -z)a{z - U2,uj2 - 2) . (2.20) 

Using (12.181) this last equation and the unitarity constraint become 

X^ -|- \ x" 

aB{z)aB{z - U2) = ',,,/, , aB{z)crB{-z) = 1. (2.21) 

x~^ + l/x~*~ 

This equation is solved in Appendix [Dl 

The above derivation was not very careful in treating the square-root branch cut in 
(I2.18p . so there may be a sign ambiguity on the right hand side of the first equation in (I2.2ip . 
The equation as written is in fact correct, to see that recall that the crossing equation was 
originally written for 5*0 in [32] as 

5*0(21 + U2, 22)5*0(21, 22) = l//(Zl, Z2) 
5*0(21, 22 - ^2)5*0(2:1, 2:2) = l//(2i, 22) 
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and the dressing factor satisfies a very similar set of equations 



X 



(2.23) 



a[Zi,Z2 -UJ2 (J{Zi,Z2 



Xi 



hf{Zl,Z2) 



where in both cases 



f{Zi,Z2) 



1 2 I ^ A ^1 



Applying the monodromies twice on fl2.2Qp gives 



(Tb{z asiz -UJ2 



So{z,-zYao{z,-zf 



(2.24) 



(2.25) 



so indeed (Jb{z) should satisfy (12.211) 

The solution found in Appendix [D] following the methods of ^0] is 



<Jb[z 



P«XB{a:+)-«XB(2-' ) 



where 



Xb[x) 



dz 1 , sinh27rqf(z + 1/z) 

log — —- 

2711 X — z 27Tg{z + l/z) 



(2.26) 



(2.27) 



A more general solution to the crossing equation is gotten by multiplying by the exponent of 
an odd function of n, so ase-^"'^'^^^^ is also a solution. In particular the simplest modification 
is a linear functioqj ase*"^. Such a linear shift does not seem to be consistent with the one- 
loop calculations in |2Tj. All that follows assumes that this is indeed the correct boundary 
dressing phase, which should be verified by further direct calculations. 

The leading behavior of as at strong coupling can be extracted from the first term in 
the expansion (JD.22P 



(^b(z 



exp 



2ig I X 



x^ 



log 



X 



X"" 



2ig I X 



X 



log 



X 



X 



At leading order at strong coupling x^ ^ 

a sip) ~ exp 



o±ip/2 



SO 



P 



Aig cos - log 



sm 



sm -■ 



Together with the usual dressing factor restricted to the boundary l/a{p, —p) 
reflection phase from the boundary at strong coupling is 



Ro{p) ^ exp 



p , p p 1 — sin f 

8ig cos - log cos - + Aig cos - log 

2 2 2 1 + sm - 



2 J 



(2.28) 

(2.29) 
the total 

(2.30) 



^Recall that x^ /x = e*^, so if one uses as as the definition of the square root in (|2.18p . then indeed 

So{z,-z) ^e'P(TB{z)/cr{z,-z). 
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The analogue calculations can be done in string theory by scattering in the classical solution 
of [21] along the lines of [28] . In addition to the usual sine-gordon scattering along the sphere 
part of the sigma model one needs to include a sinh-gordon contribution from the AdS part. 
The derivation of the boundary scattering phase above and in the appendix is for a funda- 
mental magnon. For a bound state of Q magnons it is defined as Rq{z) = a^{z)/a^'^{z, —z). 
It is evaluated by considering the Q{Q — 1) scatterings of the constituent magnons of each- 
other and the Q reflections of the constituents, a^'^ is then the product of Q^ fundamental 
dressing factors and is the same as the usual bulk bound state dressing factor, a^ is the 
product of Q fundamental boundary dressing factors which ends up being identical to fl2.26p 
with the appropriate x"^ for the bound state. 

2.5 Twisted boundary conditions 

To construct the finite size open spin-chain one has to impose boundary conditions at the left 
end of the chain as well. The left boundary conditions are also associated to a Wilson loop, 
so are very similar to the right boundary conditions jj If the Wilson loop is straight, such 
that the angles 9 = (p = 0, the boundary conditions are completely compatible, and the left 
reflection matrix is identical to the right reflection matrix[] The resulting open spin-chain 
is then described as a single p5u(2|2) spin-chain with periodic boundary conditions and a 
symmetry requirement, that for every magnon of momentum p on the original segment there 
is another magnon with momentum —p on the mirror segment, as was also the case for the 
SU{2) sector in ^. 

When the angles 6 and (p are different from zero, the two segments of Wilson loop to which 
the local operator is attached are not aligned. This means that the boundary conditions are 
not the same, but have a relative rotation with respect to each other. In the identification of 
PSu{2\2)l and psu{2\2)ji there would be different matrices instead of the a^ appearing in fl2.3p . 
This can be implemented on the reflection matrix (12.111) by a f/(l)^ C SU{2\2) rotation. For 
example choose the rotation to act on the fundamental representation of Psu(2|2)r by 

wi -^ e'%1 , W2 -^ e-'%2 , ^1 -^ e''*'^i , ^2 ^ e-'^^2 ■ (2.31) 

This is just the representation matrix of the spacial rotation by angle and i?-rotation by 
angle 9. This twist matrix is labeled G such that the reflection matrix from the left boundary 
becomes 

•")-(p) = G^§-(-p,p)Gf. (2.32) 



^That is not required, of course, and it would be fun to consider open spin-chains with boundary con- 
ditions associated to different objects in the gauge theory, describing for example local insertions at the 
endpoint of an open Wilson loop, as in |41j . 

^The phase / also matches on this wall, since the total momentum including the magnons and their 
images is zero. 
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Figure 7: A graphical representation of the boundary Bethe-Yang equation in the 
open system with two boundaries (a), or in the doubled pictures (b) with periodic 
identification. 



This is be very similar to the twisting arising in the cases of /3 and 7 deformed A/" = 4 
SYM [l2]. These models were discussed recently from the point of view of the TBA and 
the F-system in [13], HU HHl HHl HTj. The twist appears in the TBA equations as chemical 
potential terms. 

It is now possible to write the boundary Bethe-Yang equations for this system. It states 
that each magnon, if scattered with the others, reflected off the right boundary, scattered 
back all the way to the left boundary and then back to its original position (see Figure [7]), 
it would pick up a trivial phase. In matrix form for a spin-chain with L sites they are 



M 



i-1 



^2iLp^ = Y[ §(p„P,)M(^Hp.) n ^(Pi' -Pi)^^''\-P^)ll^{P^,PJ) ■ (2-33) 



j=M 



As usual this set of equations can be diaganolized by using the nested Bethe ansatz equations, 
as was done for the case of the D5-brane defects in [31]. The set of equations for the case at 
hand is almost identical, one just has to insert the extra phases from G, which acts diagonally 
on the nested equations. 

In addition to the fundamental magnons the spin-chain has bound states. In the physical 
theory they are in totally symmetric representation of each of the p5u(2|2) algebras. The 
Q^^ symmetric representation can be written in terms of homogeneous polynomials of degree 
Q in the two bosonic and two fermionic variables Wi, W2, ^3 and '^4. The representation is 
comprised of 



Bosons: 
Fermions: 



I Q-l 



/ = o,---,g, 

/ = 0,--- ,Q-1 



"^^'"'^3^4; / = 0,---,Q-2 



W^W2 



I Q-l-l Q 
W^W2 V4 



1 = 0, 



• Q-l, 

(2.34) 
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The states of the mirror model are of more importance in what follows. They are in the 
totally antisymmetric representations [H], which are homogenous polynomials in a pair of 
fermionic and a pair of bosonic variables with the opposite lables '&i, '^2, w^ and w^. The 
states and their transformation under G are 



Bosons: 



Fermions: 






«„,,Q-i-' 



(2.35) 



WnW 



3 "^4 



^^ ^ g^(2Z+i.Q)</,+i^ u;^w;«-i-'^i ; / = 0, ■ ■ ■ , g - 1 

wiwf-'-% ^ e^(2'+i-Q)<^-'? wiw^-'-% ; / = o, ■ ■ ■ , g - 1 

Then it is easy to calculate the supertrace 

Q Q-2 Q-l 

sTrg G = J2 e^''-'^)*'^ + J^ e^^'^^'"^^"^ - 2 cos ^ ^ e^^l+l-Q)^'^> 

1=0 1=0 1=0 

Q-l 

= 2(cos - cos 6) V e(2'+^-«)"^ = 2(cos - cos ^) — 
^^-^ sin( 

z=o 



(2.36) 



3 Wrapping corrections 

An elegant way to formulate the Bethe-Yang equations is in terms of the transfer matrix, 
capturing the monodromy around the spin-chain. In the case of open spin-chains the analog 
quantity is known as the double-row transfer matrix [49] defined as 



T{z\zi, ■■■ ,zm) = sTr S{z, z^) ■ ■ ■ S{z, zmW^\zP{zm, -z) ■ ■ ■ §(zi, -z)R^^^\z - 002] 

(3.1) 
where the trace is performed only over the states of the magnon with generalized rapidity z 
and it carries the matrix indices of all the other magnons.This is illustrated in Figure El 

In terms of the transfer matrix the Bethe-Yang equations take the very simple form 
T(2;i|2;i, ■ ■ ■ , zm) = — 1- 

Within the context of integrability of A/" = 4 SYM the double row transfer matrix was 
used in calculating finite size corrections to the spectrum of insertions into determinant 
operators (giant gravitons in the dual string theory) in [501 IST]. generalizing the pioneering 
work of Bajnok-Janik [33]. Here the same is done for the open spin-chain model related to 
Wilson loops. 

It is easy to show that in the case studied here, where the reflection matrix is proportional 
to the bulk scattering matrix and commutes with the twist matrix, the transfer matrix 
factorizes to the product of two twisted transfer matrices of psu(2|2) as in Figure IHl Instead of 
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Figure 8: The double-row transfer matrix evaluated on a state with two physical 
magnons of momentum pi and p2 (a). The auxiliary magnon, of momentum p 
scatters off them, reflects off the boundary, scatters again and then does a crossed 
reflection back to the original starting point, where it can be traced. In the doubled 
version (b) there are a pair of auxiliary p5u(2|2) magnons of momentum p and —p. 
They scatter off all the physical magnons and their images as well as off each other. 
Once in the forward direction and once doubly crossed. 



writing formulas with intractable indices, the reader should be convinced by Figure dUl The 
twisted transfer matrices of the closed spin-chain were studied in |15], see also |l3l HH H6| H7] . 
The main example of interest is the ground state energy, so we can omit all the magnons 
except for the auxiliary one being traced over. The result is a transfer function Tq (z) for 
an auxiliary bound state of Q magnons with rapidity z and with the twist angles explicitly 



spelled out. Using 
fl^:^ one finds 



?Wf 



GM^^''(— 2;)G and the same manipulations as in (12.151) and 



T{z) = sTr R^^\z) R^^'^'iz - 002) 
_ Ro{z)Ro{-z + UJ2) 

So{z, -z)So{-Z + UJ2,Z- UJ2) 



sTr 



,iB)' 



-z + U2) 



z,-z 



;lE^{-z,z)G^ 



(3.2) 



(TB{z)aB{-Z + U2) 



(sTrG)^ 



The final expression is very simple because the twist matrices G commute with the scattering 
matrix and thus lead to the factorization into two traces of the twist matrix and a scalar 
phase. 

The auxiliary particle can be in any representation. For the purpose of calculating the 
wrapping effects it should be a bound state in the mirror model, which are in the completely 
antisymmetric representation. Each of the supertraces is performed in the Q*^ antisymmetric 
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Figure 9: An alternative description of the double row transfer matrix, see the 
sequence of steps in Figure [TOl The result in the doubled picture (6) is a pair of 
regular psu(2|2) transfer matrices. The only remnant of the boundary reflections 
is the phase factor a^ at each boundary and possibly a twist matrix. In the open 
spin-chain picture (a) the boundary reflection is completely diagonal. 



representation of psu(2|2) which using f l2.36p leads to 



(sTrQi 



4 (cose/) — cos( 



,sm 



sm 



(3.3) 



The ground state Z'^ has classical dimension A = J. The first correction comes from the 
leading wrapping effect, which is 



^^--^E/ ciplog(l + 7;(^'^)(.+ f)e-^^'5) 



(3.4) 



Here z = Z+UJ2/2 is the generalized rapidity for a magnon of the mirror theory of momentum 
p. L is the length of the world-sheet and is equal to the number of sites J on the spin-chain. 
Finally Eg = logx^/x~ is the energy of the bound state in the mirror theory. A similar 
expression arises from the full thermodynamic Bethe ansatz treatment of the problem, where 
Eq is replaced by sq, the solution of the TBA equations, accounting for entropy, in addition 
to the energy. In that case L is interpreted as the inverse temperature. The {x~ /x'^)'^ term 
in (13. 2 p can be absorbed in a shift L — )■ L + 1 

The contribution from the boundary dressing phase is 



aB{z + ^)aB{-z + 



ijJ2 ' 



^b(2 + ^) 



piXs(a;+(2))"iXB(l/a: (z))-iXfl(l/a;+(z))+iXs(a; (2)) 



GbKZ 



CJ2 
2 ' 



(3.5) 



since a shift by uj^j'^- replaces x — t- 1/x and by —bj^jl replaces x'^ — )■ 1/x^. For real z one 
has 1x^(2)1 > and to define x^i^jx^i^zy) one should use the monodromy property (ID.9p . 
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Figure 10: The sequence of Yang-Baxter and crossing/ unit arity relations a ^ b -^ 
c —> d used to relate Figure |8] and Figure El Recall that up to a scalar phase and 
a twist with which it commutes, the reflection matrix is the same as the regular 
p5u(2|2) scattering matrix. 



this gives 



'JBiz + ^) 



MXBix+)+XB{x')). 



(27r^)2(x+ + l/x+)(x- + 1/x-) 



asiz 



Co'2 
2 ' 



sinh(27r(yf(x+ + l/x+)) sinh(27r5f(x~ + l/x~)) 
sinh (27iu) 



(3.6) 



Where the last equality was written using the relation glx"^ + l/x"^) = u ± iQ/2 and the 
periodicity of the hyperbolic sine function. Note that the dressing phase contribution has a 
pole at n = 0, which is also p = 0, representing the contribution of zero momentum states 
for all values of Q. 

Usually one expands the log function in (13. 4p to linear order and integrates over p. In the 
case at hand this approximation is invalid, since Tq (p) has the double pole at p = 0, and 
thus is not uniformly small for all p (even for small g or for large L). Instead write 



Tc 



T, 



pole 



Q _|_ rpncg 
p2 "^^Q ' 



and use 



dp\og\l + -^] = TTC 



(3.7) 



(3.^ 
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to get 



6E 



2tt 



Y^ / dplog 1 



T, 



pole 



p2 



-2LEq 



-^gy„ *i°^(r 



l + Tae 



~2LEr 



r^pole^_2LEQ/p2 



■^EA 



ipole 2LE, 



Q 



Q=l 



rpreg^-2LEc 



^y dp — -^ — - — + 



(3.9) 

The pole term dominates at large L and at weak coupling, since it comes with e~^^'^ rather 
than e~^^^'3. And it is given by the simple expression 



T-V2^^Q=2^^^t^^sing0res 



smi 



3«(XS 



(x+)+XB(.-)) 27rv/n^ + gV4 /x 



-1)*^ sinh(27r'u) \x 



L+l 



(3. 



0) 

One may be concerned about the choice of sign on the right hand side. If the integral in 
( 13. Sp is regarded as a real integral of a positive definite quantity, then one should choose the 



positive branch of -i /Tq° °e~^^^<5 . It is actually more natural to take it to be an analytic 
expression and then the hyperbolic sine in the denominator is sinh(27r(n ± iQ/2)). See a 
careful discussion of such signs in [53] . 

By using the explicit expressions at small p 



P I 16(7^ 



' 'l + ^ + 0(p3) 



x~^ 



16g2 

p2+Q2 



163^ 



equation (13.101) becomes 



rr^e-^^. = 2^"^'^-^"^^ sin 



Q 



-i)«g 



sm( 



r)2 



1 + ^-1 



+ 0{p 



^2\ 



1 + ^ + 1 




L+l 



(3.11) 



oi{XB{x + )+XB{-->: )) 



(3.12) 



At weak coupling the dressing phase starts contributing at order g"^, so it can be ignored 
when considering the first two terms 



T.ole^-2LE, = ^ COS<l>-COSe ^.^ 



■Q 



smi 



;-i)' 



^^^^ - ^(L + 2)^^^^ + 0(^^(^+=^)) 



Q 



2L+1 



Q 



2L+3 



(3.13) 
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And then from (13. 9p . the leading contribution to the ground state energy at small g is 
on. ~ \^y ) ^\^A 2-^ n2L+i 



sin0 ^^ Q2 

^'^y^'^-T' (L^-^i(-e^^) - Li....(-e-^)) (3-14) 

2% sm0 ^ ^ 

(-167r2^2)L+i cos - cos ^ f (I) 1 

-D2L+1 



47r(2L + l)! sin0 '"^' VStt 

Where i?2L+i are Bernoulli polynomials. For L = it is Bi{x) = x — |, so 

£ = ,£ = 29^22£izi^^ + 0(9*) (3.15) 

sm0 

This is the same as the one loop perturbative calculation, equation (II. 5p . as calculated 
originally in [T5] . 

For L > this calculates the leading correction to the energy of the ground state Z^ in 
the spin-chain, which happens at order L + 1 in perturbation theory and should hold up 
to order 2L + 2, where double wrapping and the second term in (13.91) start to contribute. 
For L > this would involve gauge theory graphs like in Figure HJi, which so far not been 
calculated directly. 

For large L (scaled with g) this should give the correct answer to all values of the coupling 
and match with classical string solutions similar to those in [21] . 

4 The twisted boundary TBA 

The calculation so far enabled the derivation of the one-loop result for the generalized quark- 
antiquark potential. To find the answer at all values of the coupling requires to solve exactly 
for the ground state energy of the open spin-chain with twisted boundary conditions. This 
can be done by using the boundary thermodynamic Bethe ansatz (BTBA) equations [51] . 

The idea is to exchange the space and time directions and instead of calculating the 
partition function (or Witten index) of the finite size system over an infinite time, to calculate 
the partition function of the theory on an infinite circle over a finite time [55]. In the case 
of a periodic model the result is a thermal partition function 

Z = Tre-«^^ = Tre-^«/T^ T = l/L, (4.1) 

where Hi is the hamiltonian for the original model on the interval of width L and H^ is 
the hamiltonian of the mirror model on the interval of width R. 

In the case of the open spin-chain the boundaries get replaced with initial and final 
boundary states \Ba) and one calculates the transition amplitude between boundaries of 
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Figure 11: The boundary TBA: In the original setup (a) one has to find the exact 
ground state spectrum for finite size L. Instead, in the BTBA approach (b) one 
considers the system on a very large cylinder of circumference R, where the spectrum 
can be evaluated exactly and propagates a boundary state over a finite time L. 



type a, (3 

Z^^p = Tr e-^^^("-« = (S^|e-^«/^|S„) , T = 1/L , (4.2) 

see Figure [TT] In general there are many different boundary conditions for the open spin- 
chain related to strings ending on different objects in AdS^ x §^. The ones addressed here are 
those related to insertions into Wilson loop operators as discussed in the previous sections 
and the two boundaries are identical up to the rotation by the angles and 9. 

The boundary state \B) is a superposition of states with different numbers of magnons of 
arbitrary momenta, subject to the symmetry under p -> —p. The amplitude for emitting a 
pair of magnons by the boundary is represented by the matrix K{p) which is closely related 
to the reflection matrix R. It is the charge conjugate of the reflection matrix analytically 
continued to the mirror kinematicq^ 



K 



aabbf ~ 



z) = C'C'Rl'^i 



■^~ + f), 



(4.3) 



where C is the charge conjugation matrix. 

In addition, the boundary state can have single magnon contributions subject to the 
constraint p = 0. These arise when the scattering matrix has a double pole at p = 0, which 
is indeed the case here. The amplitude for the emission of such a magnon is g°'°', related to 
/sTby p] 

ResK^^^^p) = 2ig''^g^'' (4.4) 



p-^O 



*In this section quantities are written in terms of real mirror rapidity z. 
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With these two quantities, and under the usual assumptions of integrabihty, the boundary 
state can be written as [57] 



\B)=M(\^g^'A\M\^^V 



d~zK-'^\z)A\,{-~z)A\Az) 



|0) 



(4.5) 



where ^aikv) is a creation operator for a magnon with quantum numbers a and d in 
p5u(2|2)i X psu(2|2)ij and A/" is a normahzation constant. This expression was written 
in terms of "out states" , where the momentum of the right magnon is larger than the left 
magnon. In terms of "in states" the same expression can be written by virtue of the boundary 
crossing relation fl2.14p as 



|S)=Ar 1+^-aL(0) exp 



■wi/2 



d~zK'^'^'\-~z)Al,{l)A\A-~z) 



t 



^L. 



|0) 



(4.6) 



Note that the in addition to fundamental magnons these expressions also includes bound 
states of the mirror model. If one wishes to be more explicit about that, it is possible to add 
an index Q to g, K and A and then sum over Q. 

The difference between the initial state and the final state, which can be labeled IjB*^^^) 
and \B^^'') is the same as the difference between the respective scattering matrices — it is a 
global rotation. 

Indeed, as mentioned before, when the two boundaries are of the type associated to the 
Wilson loops, then by repeated use of the Yang-Baxter equation and unitarity, it is possible 
to replace in all calculations the reflection matrix by a diagonal matrix with only the scalar 
phase and the twist. Then one gets 

where K is the conjugate of K, the amplitude of absorption of a pair of magnons into the 
final state. Since the order of the operators is reversed, in the "in states" basis the amplitude 
is K{z) and in the "out states" basis it is K{—z). 

The vacuum transfer matrix Tq is then the inner product of this two-particle contribu- 
tion 






K 



(L) 



l)K^R) 



bbaa^-^ 



\2, 



-5+f) 



= (sTrG)VB(5 + f )aB(-5 + f ) = (sTrG)2AQ(^) ' ' 

where after absorbing the factor of {x~ /x^Y i'^ (13. 2p by the shift L — > L + 1, Ag is given by 

(Elfl 

(2vr)2(n2 + gV4) 



A. 



Mxb{x+)+Xb{1/x-))\ 



sinh (27rM) 



(4.9) 



^Since mirror kinematics, where \x | < 1, are used in this section, then compared to the notations in the 
previous section one replaces x~ — > l/x~ . 
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A crucial ingredient in deriving the result in Section |3] is the fact that the scattering 
matrix has a pole at zero mirror momentum, seen here in Xq. This represents contributions 
from single zero momentum states and the factor g°-°' in equation (14 .Sp is related to the 
residue. For a generic bound state of Q magnons it is 



r, Jl + ^-1 

aa _ Qaa/ -^\Q V V V g«Xs(a;+)+«XB(l/a: ) ('4 10) 



The derivation of the boundary TBA equations follows the standard procedure [SI] of 
expressing the boundary state as a sum over eigenstate of the mirror hamiltonian, introducing 
densities for particles and holes and minimizing the free energy, accounting for the entropy 
of the states. As explained, the boundary state is completely symmetric in the exchange of 
p — )■ —p and under the replacement of psu(2|2)/j and psu(2|2)/,. The density of momentum 
carrying particles can therefore be defined for positive p only, representing a pair of particles 
(or holes) at momentum p and —p. At the nested level the replacement of n — ?■ —u is 
accompanied by replacement of the right and left groups. 

The mirror theory is exactly the same as in the case of periodic spin-chains and the 
Bethe equations (which lead to the kernels in the TBA equations) are also the same, except 
that one has to account for the symmetry of the states overlapping with the boundary state. 

The resulting BTBA equations are therefore the same as for the closed spin-chains [581 

60] with the following modifications: 

1. The momentum carrying densities eg are defined only for positive u. 

2. Since now the eg densities represent pairs of particles of opposite momentum, all ker- 
nels coupling to them are doubled K'^'^ {u,u') — )■ K^^ {u,u') + K^^ {u,—u'). The 
integration domain for u' is [0, oo). 

3. All other densities are symmetric under 

it\u) = el~Ji-u) , eil]^iu) = eill,i-u) , e%iu) = el"l{-u) . (4.11) 

4. The twisting matrix G introduces chemical potentials, as in the case of the usual close 
spin-chain twisted TBA. See [l3llMli5]. 

5. There is an extra driving term from the boundary dressing phase contributing to the 
equation for eg. This is equal to — logAg. 

6. The ground state energy is 

*^ = -sE/ *£log(lH-A«e-) (4.12) 
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An alternative formulation extends the definition of eg to negative u by eQ{—u) = €q{u) 
and then integrate u over (—00, 00) with the regular kernels. In this formulation the BTBA 
equations are identical to the usual twisted periodic TBA equations except for the extra 
driving term — logAg for Yq. In the notations of [60] the equations are 



logF^"^ 



r{a) 



-2^^ + log (1 + l/i;i|i) ^ K^M + log ^^1^ ^ K 



log YZ\m = -2# + log (1 + l/nS^) ^ ^<NM + log / 7,) ^ i^M - log (1 + Fg) ^ i^«^ 



+ 1A + 
1 + i/ri") 



IokK 



(a) 



iog(i + yQ)*A-r'+iog 



QM 



1 + 1/f/"^ 



vw\M 



1 + 1/F^^'^^ 



^is: 



M 



(4.13) 



log yg = -2tQ<f) - (L + 1)Eq - log Xq + log (1 + Yq,) * /C^ 



Q'Q 

m 



a=± 



log (1 + I/^Sm) * ^*^? + E l°g (1 + ^l^±^) ^ ^^' 






with 



Yn 



Y 



(a) 



(a) 



vw\M 



p vuj\AI 



Y 



(a) 



w\M 



(a) 
(D m\M 






M) 



e y 



(4.14) 



The definitions of the convolutions and all the kernels are given in [60] (though note the 
slightly different conventions from this manuscript). 
The energy is then given by 



E 



00 

-y 

47r f-^ 



du ^ log(l + Yq 
du 



(4.15) 



One can derive the simplified (and hybrid) TBA equations as usual by acting with the 
inverse kernel [K + 1)^\/ (see [HD]), which removes the \q terms from the equationsjl^ due 
to the relation \q_i{u)\q^i{u) = \q{u — i/2)XQ{u + z/2). The derivation of the F-system 
equations [61] then follows as usual [60], and they are identical to the regular equations for 
the spectral problem in AdS /CF1. 

Assuming Yg = one finds the asymptotic solution for the auxiliary particles which is 
identical to the case of twisted periodic TBA [17] 



F 



(±)o 



■w\M 



sin M6' sin(M + 2)6^ 
sin^^ 

(±)o _ COS0 



F 



COS 6' 



Y 
Y 



(±)o 



sinM0sin(M + 2)0 



vw\M 
(±)o 



sm 



coscp 
cos6' 



(4.16) 
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Except possibly for the usual subtleties in the equation for Yi. 
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Feeding this back into the equation for Yq and including the boundary driving term gives 

YS = 4(cos</. - cos^)2^^^^ AQe-(^+^)^« (4.17) 

sin (j) 

Note that the assumption of small Yq is true only for u ^ and this is not a good approximate 
solution at M = 0. Still as shown in the previous section, the pole contribution does give the 
correct one loop gauge theory result. 



5 Discussion 

This paper presents the first step in finding the exact quark-ant iquark potential in A/" = 4 
SYM: The set of twisted boundary TBA equations whose solution is conjectured to give the 
desired function. As explained in the introduction, with the extra two parameters 6 and 
the generalized potential is the same as the generalized cusp anomalous dimension. 

The resulting quantity gives the conformal dimension of cusps in arbitrary Maldacena- 
Wilson loops in this theory. The discussion here has been on the most symmetric cases, of 
the antiparallel lines or infinite cusp. But the divergences are a UV property and depend 
only on the local structure of the Wilson loop. Knowledge of the exact generalized cusp 
anomalous dimension allows therefore to renormalize any Wilson loop (with the usual scalar 
coupling) in this theory with arbitrary number of cusps. One follows usual normalization, 
now that the UV behavior near the cusp is under control, see [TT] . 

The discussion here was for real angle as appropriate for the Euclidean theory (or a 
Euclidean cusp in the Lorentzian theory). There is no reason to restrict to that. All the 
expression are analytic in and a Wick-rotation (p = iip gives the cusp anomalous dimension 
for an arbitrary boost angle yj. At large ip the result should be proportional to c/), with the 
coefficient being a quarter of the universal cusp anomalous dimension 7cusp/4. In particular 
for large ip the BES equation for the cusp anomalous dimension [62] (and its solution [63]) 
should be recovered from the twisted boundary TBA equation. 

For large but finite y? this gives the regularized cusp anomalous dimension which plays 
a role in scattering amplitudes in A/" = 4 SYM when imposing a Higgs-VEV regulator as in 

[63 EH [7]. 

The TBTBA equations contain much more information than that, though. As usual, 
by contour deformation one can calculate the exact spectrum of excited states of the rel- 
evant string (or cusped Wilson loops with local operator insertions). For large L one can 
approximate the solution by solving the boundary Bethe-Yang equations instead. 

In the case of large imaginary angle, these are excitation of the string with lightlike cusp 
describing scattering amplitudes in AdS. This spectrum is a crucial ingredient in the OPE 
approach to lightlike Wilson loops as discussed in [661 EZ] and should be closely related to 
the excitation spectrum of the spinning string as studied in [SS] (see also 
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There are many other questions left for the future. One is numerical solutions of the 
TBTBA equations, giving the interpolating function for the quark-ant iquark potential and 
other interesting quantities in this model. 

Another is to try to solve these equations at large coupling and reproduce semiclassical 
string theory results, as was done from the F-system in [7D] . Normally these techniques work 
only for large L and all other charges also large, since otherwise there is no semiclassical 
string description. In this case, though, the classical string description exists already for 
L = as the usual description of Wilson loops in AdS^ x S^. It should therefore be that the 
algebraic curves describing the string duals of these Wilson loops [ISl [HI [151 [iHl [6] can be 
derived from these equations. 

The same tools used here can be used to study other open spin-chain models which arise 
in the AdS /CFT correspondence. The simplest are the D5-brane defects, with a very similar 
symmetry to the Wilson loop. There are two natural choices for the vacuum there, one with 
and one without boundary degrees of freedom. The latter should be essentially the same as 
the model studied here, while also the former should not be much different. One would need 
to calculate the boundary scalar factor and then always carry through the two boundary 
excitations, so similar to calculating excited states in this model. 

Another interesting system to study is ABJM theory |7T]. The 1/2 EPS Wilson loop is 
known in that theory [72] , but the quantity calculated here — the quark-antiquark potential 
or cusp anomalous dimension has not been calculated in the gauge theory. The leading 
classical result from string theory on AdS^^ x CP^ is the same as for AdS^ x S^ (with the 
appropriate identification between the string tension and the gauge coupling). 
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A Supersymmetry 

The vacuum of A/" = 4 SYM is invariant under the P5'f/(2,2|4) superconformal group, 
summarized here following [55] . 
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Denote by -L"g, X"g the generators of the SU{2)l x SU{2)ji Lorentz group, and by R ^ 
the 15 generators of the i?-syninietry group SU{A). The remaining bosonic generators are 
the translations P^'^, the special conformal transformations K^a and the dilatation D. 
Finally the 32 fermionic generators are the Poincare supersymmetries Q^, Q and the 
super conformal supersymmetries S^, SaA- 

The commutators of any generator G with X^g, X"^ and R ^ are canonically dictated 
by the index structure 

[L-^, G^] = 5^G" - |<5^G^ , [L'^^, G,] = -6^G^ + i^^G, , 

[£"^, G^] = S^G'' - Id^G"^ , [L°„ G,] = -6-G, + i5|G, , (A.l) 

[R^^, G^ = 5gG^ - i^^G^ , [i?^^, Gc] = -6^Gb + i^^Gc • 

while commutators with the dilatation operator D are given by [D, G] = dim(G) G, where 
dim(G) is the dimension of the generator G. 
The remaining non-trivial commutators are 



{Q%Q''}=6iP^\ 




{5q,, 5^5} — 


^iJ^a/3 ' 


.-^a/3' Qa\ = ^Z^pA ' 




[i^«,,Q^1 = 


= ^^^^ , 


[P'^^ S^] = -b^^Q^^ 


5 


p°^s^A] = 


^ --^-Ql , 


{Q%S^,] = blL-,-.b-,R\ 


+ i^f 5|^ , 




{Q-\s,^}=biL%- 


-blR 


^ + i5|5°D , 





(A.2) 



To write down the integrable spin-chain for A^ = 4 SYM one chooses a vacuum corre- 
sponding to one complex scalar field, usually labeled Z = $^-|-i$^. This choice of spin-chain 
vacuum breaks the symmetry group PSU{2,2\4:) -^ PSU{2\2)'^ with supercharges (Qq,5'") 
and (Qf , "S"^) where a,b E {1,2} and d, 6 G {3,4}. Explicitly, the two copies of psu(2|2) are 
the following subset of the generators of p5u(2, 2|4) 

i^a f^a QO. Qa ja ja pa r>a 1 A*^ R"^ • 

^a—^a^ '-'a — '-'a^ ^ (S — -'-' (3 ^ ^ b — ■'^ b 2"b-"'c) 

_.; _ ._ _ _ _ A.3 

In the presence of the Wilson loop or domain wall this is further broken down to a single 
copyof P^?7(2|2). 

B The p5u(2|2) spin-chain 

The spin-chain description of single trace local operators, which we will use also for the 
purposes of studying the boundary changing operators involves choosing a ground state, 
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usually taken to be Tr Z'^ with Z = $^ + i$^ and considering the excitations about it. The 
Choice of scalar Z breaks the symmetry group PSU{2, 2|4) — )■ PSU{2\2)'^ x f/(l). Magnons 
are therefore classified by representations of the broken group. The basic magnons are in 
the fundamental representation of each of the PSU {2\2), so we can first treat the magnons 
as if charged only under one of the groups, remembering to pair them up later. 

When constructing magnon excitations on the spin-chain it is useful to consider the 
central extension of the p5u(2|2) algebras. The commutators of L'^a, L"-, i?"^ and i?". are 
inherited from (lA.ll) . The central extension for the Q and 5* commutation relations are 

{Q:, gf } = e"^e„,P , {SI 5j} = e^pe'^'K , 

and likewise for the second copy of psu(2|2). 

The fundamental representation has a pair of bosons </>'* and a pair of fermions ^". The 
algebra acts on them by 

jx 

where a, the spectral parameters x"^, the coupling g and the magnon momentum p are all 
related by 

x+ + — - - X" = - , a^ = ig(x~ - x+) , e*^ = — . (B.3) 

x+ x^ g x~ 

The central charges for this representation are given by 

C = a2-i, P = gfil-e'n, K = g r\l - e^'^) . (B.4) 

For a single magnon we can eliminate the parameter / by rescaling Q — )■ /^'"^Q, S — t- 
j-i/2^ and ijj — )■ f^''^ip. It is crucial though for integrability when constructing muli-magnon 
states. In that case we take for the k^^ magnon /^ = e^^j=k+iP3 ^ Jt is simple to show by 
induction that then for M magnons the total central charges are 



M M 

■ 
fc=l fc=l 



P = y^P, = gfM{l - e^S-i^'^) , K = y2Kk = gfllil - e-^S"-^'^) . (B.5) 



The S-matrix exchanges two fundamental representations and has the general form 



2 

31 



Up to an overall scalar phase, 5*0, all the terms in this matrix are fixed by symmetry |73[[71]. 
This is achieved by imposing that the left and right hand side transform in the same way 
under psu(2|2). The solution is 

12 _ c -^2 ~ ^1 ol2 _ /ll2 I 1 _ 2 /"^^ ^1 ^2 ~ ^1 

rf \ I I I fY>' ry ry'^ ry 



A'' = So i i , B'' = A''{l-2 

D'^ = -So , E'' = D 



12 c 1712 nl2 / 1 <-, J^/"''2'^l -^2 -^1 



\ I rr> nr' /y rr>' 

I 1 1 2 1 

^12 _ c -^2 ~ ^1 ^12 _ c ^0-10-2 J2 1 ^2 ~ ^1 m ^^ 



X2 - x^ 12 c 2«ia2 1 1 x^ 



i^ — 'JO^ + ) -T — oo 



+ ' ^^^+^+^-f. 1 1/^-^-^- ^+ ' 



xfx^ gf2l - ^jx^x^ X2 



Xo Xi Ju'\X'-)yJ.I'2-^ X/J/qXiXo X 



i7^^ = Ki2 = 5o^^i^^ 



5' ^2 -^1 

C Reflection matrix 

The refiection matrix is M(p) = Ro{p)E>{p,—p), where S is the scattering matrix written 
above with the factor of So removed. Ro{p) = o'b{p)/(t{p, —p) is the boundary scalar factor 
discussed in Section 12.41 and evaluated in the next appendix. Under p — )■ —p the spectral 
parameters transform as x^ — )■ — x^. Therefore using the same ansatz as in (IB. 61) . the 
coefficients of the reflection matrix are 

. D ^" Rl2 p X-{{x+f + X-~) 

^12 = -Kq— T , -D = — ito" 



X'^ [X'^Y{1 + X+X" 

12 12 ^^^ ^ (r~^^^ 

D = — i?o ) E = Ro' 



x~(l + x+x~) 

rr>\' I rp ( rp'T' I /y 1 ( /y\' 'f'^ 

j^l2 _ rl2 _ D X -|- X 12 _ -D l-^ ~r -^ A-^ -^ 

2x+ ' x~(l + x+x~) 



(C.l^ 



12 r^l2 D 17I2 „-D -^ V"'' "tX JI,X X 



2x+ ' (x+)2(l + x+x-) 

D Boundary dressing factor 

D.l derivation 

The boundary dressing phase ^^(m) should satisfy the crossing and unitarity equations (12.2 ip 

x~ ~l~ 1 1 x^ 
aB{u)aB{u) = , , , , , aBiu)aBi-u) = 1 . (D.l) 

x+ + l/x+ 
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This equation can be solved the same way the bulk one is |10] (which is significantly simpler 
than the way it was originally found [751 ES]), see also [76]. Take the ansatz 

asiu) = e^(^s(^+)->^s(^")) , (D.2) 

and defining 5"b(x) = e'^^^'^^^^~^^'^^^^^^\ then using that under crossing transformations x^ — )■ 
1/x^ the crossing equation (ID .11) can be written as 

aB{u)aB{u) = aB{x+)aB{x~) = ^ ^ = ^ / . (D.3) 

x^ + l/x+ u + 1/2 

In terms of the shift operators D"^ = e^^^^ this equation becomes 

aij(x)^^-^" = {x + l/xf--^^ = n^""^^ . (D.4) 

So 

aBix) = ix + l/x)^[^] , F[D] ~ §^^ (D.5) 

Taking -F[-D] = — 1 will not lead to the desired answer, rather will introduce the false analytic 
structure. Instead one can try 

n-2 n+2 °° 

F[D] = -^ + -^ = Y: {D~'- + D-'-) (D.6) 



Thereforqi^ 



l-D-2 l-D+2 

n=l 



~ / \ TT/ • \/ • \ smhvTM ,^ , 

(Tb[u) = \\\u + in)[u — m) = (D.T) 

n=l 



This expression does not quite work. The shifts by ±| gives 

^ ^ a-(M-V2) sinh7r(M-V2)M + V2 ^ ^ ^ ^ ' 

This differs from ( ]D.4|) by a sign, which can be fixed by taking 



~ / X sinh 2Tru 

(^b{u) = -— (D.9) 

ZTTM 



This corresponds to the choice 



n- r)+ °° 

F[D] = -—— + --— = Y: {D- + D^-) (D.IO) 



n=l 



Of course one could replace the 2 in ( JD.Qf) with any even integer, but this would introduce 
extra poles. 



^^Thc product is divergent, but this is the natural ragularization of it. 
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Under crossing transformation x — )■ 1 /a; so we can interpret this last equation as the 
discontinuity of xb across the cut in the u plane between ±2g. Therefore 

f^'^^"' dw x(u)-l/x(u) 1 1 sinh27rw ,^ , 

Xb{u)= o ■ t \ TTT^ -log^ • (D.ll) 

J~2g+i0 ^^* X[W) — 1/X[W) U — W I 2tTW 

Switching to the x coordinate and ignoring an irrelevant constant gives 

. f dz 1 smh2TTg{z + l/z) m1o^ 

The function e"*-*^*^^'-"^^ has a cut for u E [—2g, 2g] in the sheet where \x{u)\ > 1. Crossing 
the cut gives |a;(M)| < 1 where 

^-ixixju)) _ ^ixa/^(n)) ^^^ (D.13) 

sinh 27rM 

On this sheet, in addition to the cut, this function has poles for all half integer imaginary u, 
except for u = 0. 

D.2 Expansions 

As is done with the bulk dressing factor, we can expand the boundary one for large x aa^i 

^XB{x) = f:'-^ (D.14) 

r=l 

Expanding (ID.12P at large x are using z = e^^ we find the integral expression 

f^'^ dtp 2ir, sinh 47r5f cos ^ /t^itn 

Cr{g)= — e^-'^log— ^— -^ (D.15) 

Jo 27r Airgcosip 

This integral can also be written as 

c.(,) = 2(-ir^y^~^|^rft (D.16) 

One can then expand at weak coupling 

oo 

c.(^) = 5^cW/"+'^ (D.17) 

n.=0 

and the explicit factors are 

"^^ n + r n!(n + 2r)!^^^+ ^ ^ ^ ^^ 



^^Onc can readily check that only even powers will appear in the expansion. 
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There is also an asymptotic strong coupling expansion 

oo 



(D.19) 



n=0 



Apart for the linear and constant terms there are only odd inverse powers of g and the 
explicit factors are 



4°^ 



\r+l 



4r2- 1 

(-ir 


2r 
16 



(D.20) 



T{n + r-\)V{n-r-\)C{2n) n = l,2, 



(47r)2"+i 
42n+i) = o n=l,2,--- 

One can now resum the series and express xb as 

oo 

ixB{x) = Y,d^''\^)g'"' 

where 

rf(o) (x) = -4 + 4 ( X + - j arccot x 



(D.21) 



n=0 



rf(2")| 



16 



X 



T{n + i)r(n - |)C(2n) 2Fi{n + i 1, | - n, -^ 



(47r)2n+i 
The hypergeometric function is a rational function of x^. 
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